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q ' Abstract. We consider the recent RHIC data on the transverse single spin asymmetry (SSA) A^, 

, measured in p'p — > 7Z°X processes at mid-rapidity by the PHENIX Collaboration. We analyze this 

f**^ ■ experimental information within a hard scattering approach based on a generalized QCD factoriza- 

tion scheme, with unintegrated, transverse momentum dependent (TMD), parton distribution and 
fragmentation functions. In this kinematical region, only the gluon Sivers effect could give a large 
contribution to An; its vanishing value is thus used to give approximate upper limits on the gluon 
Sivers function (GSF). Additional constraints from the Burkardt sum rule for the Sivers distributions 
are also discussed. 
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^ : FORMAL APPROACH 
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Transverse single spin asymmetries can originate, even with a short distance helicity 
P-i! conserving pQCD dynamics, from spin-A:^ correlations in the soft components of the 

hadronic process A^B — ► CX. According to the hard scattering approach to hadronic 
interactions developed in Refs. |Q1 0, [3D, based on the assumption of a generalized 
QCD factorization scheme which involves unintegrated TMD parton distribution and 
fragmentation functions, the general structure of the cross section for the polarized 
^ ■ hadronic process (A, Sa) + (B,Sb) — > C + X reads (see Ref. yfl for more details) 

Ecdo (A,S A ) + (B,S B ^C + X ^ f dXadxbdz 



X J{k LC ) p a ^ A f a/A j A (* a ,*La) Pxjx? 8 fb/B,S B i x b^± b ) 

x ^^^^S^+^+fl)^^^), CD 

where all parton intrinsic motions are fully taken into account, both in the soft, non 
perturbative components and in the hard, pQCD interactions. 

The main features of Eq. ((U) are the appearance of several spin and k± dependent 
distribution and fragmentation functions and the non-collinear partonic configuration 
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which lead to many k± dependent phases. In Ref. |3|] it was explicitly shown that the only 
sizeable contributions to the transverse single spin asymmetry A#(p' p — > TtX), in the 
kinematical region of large positive xp come from the Sivers [4] and, less importantly, 
from the Collins [5] mechanisms. Moreover, while the quark contribution is totally 
dominant at large, positive xp values (for polarized protons moving along the positive Z- 
axis), the gluon contribution may be sizeable in the mid-rapidity and negatives regions. 

Data in the mid-rapidity region are available from the E704 [g] and PHENIX [7] 
experiments; the region of negative values of xp has been covered by the STAR [|80 
and BRAHMS experiments Q9Q. In these kinematical regions A^(p^ p — > nX) is largely 
dominated by the Sivers effect [0], all other contributions being almost vanishing, and 
Eq. ® gives [H: 

E K do^ E^do 1 ^ f dx a dx b dz , 2 , ,2, ,3, sn m,/, n 

-^-^^LJ ^ d *±.*(*±.-fc)'(*±«) 

do ab ^ cd 

X A fa/pl ( X a, k±a) fb/p( x b,k± b ) S (x a , X b , S, t, u)8(s + t + u) t> K / c {z, k± K ) ,(2) 

where 

A fa/pl (x a ,k_La) = f a / P T [Xa,k±a) ~ f a /pl (*a, k_L a ) = ^fa/pl ( x a, k_La) COS (j) a . (3) 

A N f a / p i(x a ,k± a ) [or fyp(x ai k± a )] is referred to as the Sivers distribution function of 
parton a inside a transversely polarized (along the Y-axis) proton (moving along the 
Z-axis). a is the azimuthal angle of the intrinsic transverse momentum k± a of parton a. 

The azimuthal phase factor cos0 a appearing in the numerator of A#, Eqs. © and 
©, plays a crucial role and deserves a comment. The only other term depending on ty a 
in Eq. (O is the partonic cross section, in particular via the corresponding Mandelstam 
variable i . Therefore while at large positive xp (/-channel dominated) the integration over 
ty a does not necessarily suppress Am, for negative values of xp (w-channel dominated) 
one is roughly left with the d 2 k± a cos (j) a integration alone, which cancels the potentially 
large Sivers contribution. As a consequence, one cannot get significant information on 
the gluon Sivers distribution from the recent STAR and BRAHMS data at negative 
values of xp. (Notice that at lower values of y/s the suppression induced by the cos0 a 
dependence would be less drastic.) 

The same arguments do not apply to inclusive hadronic processes at mid-rapidity and 
moderately large pj values, for which data from PHENIX |7|] are already available, for 
neutral pions and charged hadron production. For these processes, the gluon contribution 
is dominant and the Sivers effect can survive the phase integration. The possibility of 
accessing the gluon Sivers function has been also investigated in Refs. IIo[ Til 121. 



Indirect constraints on the GSF could also be obtained from a sum rule for the Sivers 



distribution recently derived by Burkardt 111311 . The Burkardt Sum Rule (BSR) states 



that the total (integrated over x and k±) transverse momentum of all partons (quarks, 
antiquarks and gluons) in a transversely polarized proton must be zero, 

= £ (k ± ) a = fdxf d 2 k ± k ± £ Af a/p1 (x, = . (4) 



In the following we shall simply check whether or not the proposed parameterizations 
of the Sivers functions are compatible with the BSR. 

PHENOMENOLOGY 

We consider the PHENIX data [7] on A N for the p ^ p — > 7t°X process at RHIC, at = 
200 GeV, with pj ranging from 1.0 to 5.0 GeV/c and mid-rapidity values, \r\ \ < 0.35. 
In this kinematical regime, at the lowest pj values, x™ m can be as small as 0.005. 
Therefore, partonic channels involving a gluon in the transversely polarized initial 
proton dominate over those involving a quark. This gluon dominance, together with 
the (almost) vanishing of all possible contributions to Apj other than the Sivers effect, 
allows to put upper bounds on the GSF. Due to possible mixing with quark initiated 
contributions, the same is not true for the E704 data [6] at lower energies and comparable 
rapidity and pj ranges. 

In Ref. [1] we have shown that reasonable fits to the SSA for the p' p — > nX process 
at large positive xf can be obtained by using valence-like Sivers functions for u and 
d quarks, which turn out to have opposite signs. The use of valence-like u and d Sivers 
functions alone predicts an almost vanishing SSA in the mid-rapidity region, compatible 
with available data (see Fig. l(left), for the PHENIX results). Their parameterizations [Ql 
[Till are also compatible with the BSR, Eq. ©. Although a large gluon Sivers function 
would not modify the analysis of the SIDIS, E704 and STAR data at large positive xf, it 
would strongly affect the description of the mid-rapidity PHENIX data. 

In what follows (see Ref. ||15|l for details) we therefore try to understand what is 
the maximum value of | A N f g / p ^ (x,k±)\/ 2f g / p (x,k±) allowed by the PHENIX data; our 
results are summarized in Fig. 1 (left panel for the SSA and right panel for the GSF): 

• The thin, solid line in Fig. l(left), results from computing Am using only the valence- 
like u and d Sivers functions of Ref. yj]. 

• The dot-dashed curve in Fig. l(left) has been obtained by saturating (in magnitude) 
the GSF to the natural positivity bound [see Eq. ©] 

A N f g/pr (x,k 1 _) = -2f g/p (x,k 1 ). (5) 

The sea-quark Sivers functions are again assumed to vanish. This choice leads to a SSA 
definitely in contradiction with data and to a strong violation of the BSR. 

• The thick, solid curve in Fig. l(left) has been obtained still assuming that there is no 
sea-quark Sivers contribution, and looking for a parameterization of A f g / p ] yielding 
values of A^ falling, approximately, within one-sigma deviation below the lowest pr 
data. The corresponding x-dependent part of the GSF, normalized to its positivity bound, 
\A N f g ip\ (x)\/2f g / p (x), is shown as the solid curve in Fig. l(right). It leads, within the x 
range covered by the data, to a strong violation of the BSR. 

• We finally consider the inclusion of all sea-quarks (u s , it, d s , d, s, s) contributions by 
using a non-vanishing positive Sivers function which saturates the positivity bound [that 
is, A f q jp\ (x,k±) = 2fq s / p (x,k±)]. These contributions could then cancel the negative 
contribution to An of a possibly large GSF. We then look for the largest negative GSF 
which, together with a positive maximized sea-quark contribution, leads again to the 




SSA represented by the thick, solid line of Fig. l(left). This curve results now as the 
sum of the (maximized) sea and valence quark contribution (dotted curve) and that of 
the new GSF (dashed curve), which is plotted as the dashed curve in Fig. 1 (right). It is the 
largest (overmaximized) gluon Sivers function compatible with PHENIX data. Within 
the x range covered by the data, the (over)maximized sea-quark Sivers distributions give 
a positive contribution which strongly suppresses the negative contribution of the GSF, 
so that in this scenario the BSR is satisfied within a 10% level of accuracy. 

Summarizing, our analysis shows that the PHENIX data on Apj(p^ p — ► n°X) allow 
to put significant quantitative bounds on the magnitude of the GSF. Similar conclusions 
have been recently reached by studying the Sivers effect in SIDIS off a deuteron tar- 
get M- 
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